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Characteristic structures such as canals and ridges — intermode traces — emerge in the spacetime 
representation of the probability distribution of a particle in a one-dimensional box. We show that 
the corresponding propagator already contains these structures. We relate their visibility to the 
factorization property of the initial wave packet. 



I. INTRODUCTION 

The Zeeman or Stark effect in a hydrogen atom Q|, the phase change at a conical intersection in a polyatomic 
molecule and the nuclear shell structure are three striking phenomena illustrating the importance of degeneracies 
in quantum mechanics. But how can degeneracies play a key role in such a well studied problem as the particle in the 
box? In the present paper we show that the characteristic structures @-^] in the spacetime representation of a particle 
in a box shown in Figs. |l| and ^| and referred to as quantum carpet are a consequence of a manifold degeneracy [jH^. 

Our analysis relies on three observations: 

(i) The wave function of the particle at time t is the integral of the initial wave function and the well known 
propagator of the box problem jl| . 

(ii) The probability of finding the particle is the absolute value squared of the wave function. In the propagator 
formulation of quantum mechanics this translates into the two-dimensional integral of the product of the initial 
wave function and its complex conjugate multiplied by a kernel. This kernel is the product of the Green's 
function and its complex conjugate evaluated at the two different initial conditions. 

(Hi) The quadratic dependence of the energy on the quantum number factorizes in the product of the Green's 
functions into the product of two linear dependences. 

These observations allow us to rewrite the kernel in a way which brings out most clearly that the carpet structure is 
inherent to the propagator. Three results stand out: 

(i) The kernel is only non- vanishing along straight lines in spacetime. 

(ii) The steepness of these lines is discrete. 

(Hi) These lines either start at the difference or the sum of the integration variables. 

But where is the degeneracy? The first degeneracy has its roots in the dispersion relation connecting in a quadratic 
way kinetic energy and wave number. This forces many quantum numbers to contribute to the same spacetime line. 
Another degeneracy arises when the product of the initial wave function times its complex conjugate factorizes into 
the product of two functions which now contain the sum and the difference of the integration variables, only. This 
degeneracy in the initial wave function either enhances or suppresses the intermode traces resulting from the spectrum 
degeneracy. We note that the phenomenon of fractional revivals j|,|ll],|l2| is closely related to the intermode traces. 
Our paper is organized as follows. In Sec. O we cast the probability to find the particle at time t at the position 



x in the box into a kernel formulation. Within Sec. Ill we rewrite this kernel to bring out the intermode traces and 
discuss in Sec. IV the consequences of the factorization property for the emergence of the carpet. We illustrate these 
results in Sec. |V| for the example of a Gaussian wave packet. We conclude in Sec. VI by summarizing our main results 
and giving an outlook. 
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II. THE PARTICLE IN THE BOX: A BRIEF SUMMARY 



In this section we briefly summarize the essential ingredients of the problem of the particle in the box. In particular, 
we express the probability W(x, t) to find the particle at time t at position x as an integral of the product of the 
initial wave packet, its complex conjugate and a kernel. This kernel is the product of the two Green's functions at 
two different initial positions. 

We consider the motion of a particle of mass M caught between two walls separated by a length L. The wave 
function %j)(x,t) of this particle at time t reads 



ip(x,t) = ^ u m (x) exp (-iE m t/h) . 

m— 1 

Here the expansion coefficients 

tp m = / dx' u m (x') ip(x') 
Jo 

follow from the wave function ip(x, i = 0) = ip(x) at time t = and from the energy eigenfunctions 



u m (x) = Afsm(k m x) 
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with eigenvalue 
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Here we have introduced the revival time T = 4M L 2 / (irh) . 

Note that the form of the energy eigenfunctions u m together with the wave numbers k m guarantee that the wave 
function ip satisfies the boundary conditions 



ip(x =0,t) = i/j(x = L,t) = 

enforced by the two walls at x = and x = L at all times. 

The constant M = (2/L) 1 / 2 results from the normalization condition 



dx u m {x) 2 =Af 2 I dx - [1 - cos(2A; m a;)] = N 2 
a Jo 2 2 
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For further calculations it is more convenient to have the summation in the representation Eq. (Q) of the wave 
function run from minus infinity to plus infinity rather than from unity to plus infinity. For this purpose we first 
decompose the energy eigenfunctions 



u m {x) = ^\ ( 
into right and left running waves which yields 
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When we define the expansion coefficients ip m for negative values of m by 

1p-\ m \ = -1p\ m \ 

we find the compact representation 
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of the wave function. Here we have used the fact that ipo = since uq = 0. 
When we now substitute the explicit form 



1 



W2L Jo 



dx' 



exp mm — — exp —mm — 



of the expansion coefficients into the expression Eq. ([L0|) of the wave function we arrive at 

i-L 



ip(x,t) 

Here we have introduced the Green's function 
G(x,t\x') = — ]T 



dx'tp(x')G(x,t\x'). 
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of the particle in the box. 

With the help of this result we can represent the probability 



W(x,t) = \iP(x,t)f 



dx' / dx" ip* (x) ip(x") K(x,t\ x' 7 x") 



(11) 
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to find the particle at time t at the position x in terms of the initial wave function if, its complex conjugate if* and 
the kernel 



K(x,t\x',x") = G*(x,t\x') ■ G(x,t\ x") 



(15) 



consisting of the product of the two Green's functions at the two different initial positions. 

In the next section we show that the kernel already contains the intermode traces in the carpet. They emerge in a 
striking way provided the initial wave packet satisfies an appropriate factorization condition. 



III. A NEW REPRESENTATION OF THE KERNEL 



In this section we first find a new representation of the kernel. We then discuss its properties. For the details of 
the calculations we refer to the appendix. 



We first calculate the kernel defined in Eq. (I15f). For this purpose we substitute the expression Eq. (13) for the 



Green's function G and an analogous expression for G* into the definition, Eq. (|15|), of the kernel. After minor algebra 
we arrive at 



K{x,t\x',x") 
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When we recall the relations 



x 1 + x" 

m'x' - m"x" = (to' - to") h (to' + to") 



(17) 



and 



x' + x" 

m'x' + m"x" = (m' + to") h (to' - to") 



(18) 



we can cast the kernel into the form 
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K(x,t\x',x") = 



4L 2 



V 

- V 



2L ' 2L 
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It consists of four contributions, each of which involves the function 

OO s 

V(j]\ () = E exp [iw(m' + m")rj\ cxp < —i-n{m' — m") 
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Note that the four terms differ in the arguments of V. 

The function V looks rather complicated. However in the appendix we derive the equivalent representation 

OO 

v(v,0= £ (-I)"'' e*™' *[*„,,(*,*) + C] 
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where 



Xn,l(*>*) = | ~ n Y^ -I- 



(21) 
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With the help of this relation the kernel K for the probability W to find the particle at the time t at position x reads 



K( 
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Xn,l(x,t) 
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This is the key result of the present paper. We now discuss the properties of the kernel using this representation. 
Here we focus in particular on the behavior of K in the spacetime strip defined by the left and the right walls at x = 
and x = L. 

We note that the kernel K has quite a characteristic form. Due to the appearance of the <5-functions it is different 
from zero only along the four sets of straight spacetime lines 



and 



t n ,i(x;x' + x") _ 1 fx x' + x" 
T/2 ~ n\L T 2L 



t n ,i{x; x' — x") _ 1 / x x' — x" 
T~/2 = n \L T 2L 



(24) 
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The steepness nT 1 of these world lines is quantized in terms of the integer n, which due to the summation can take 
on any value from minus infinity to plus infinity. 

These lines enter or leave the spacetime strip through the left wall at x = at the times 



t n<l ( x = 0;x' + x") _ Ix' + x" I 
T/2 ~ T n 2L n 



(26) 



and 
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t nl (x = 0;x' -x") _ Ix'-x" I 

~ Tj2 ~ T n^2I~~n- (?) 

We note that apart from the integers n and I the integration variables x 1 and x" determine these times. Hence for 
one set of n and I values these spacetime lines are determined by the integration variables x' and x". However, there 
is a degeneracy in these lines: the integration variables x 1 and x" enter these expressions either in their sum or their 
difference. Hence it is left to the initial wave packet to either enforce or suppress this degeneracy. This either enhances 
or reduces the visibility of the intermode traces, as we show in the next section. 

We emphasize that there is already one degeneracy guaranteed by the quadratic energy spectrum Eq. (|]) of the 
particle in the box: In the function T> the summation indices m! and m" enter the expression in the square brackets 
multiplying time in their sums m! + m" only. Indeed by multiplying it with the prefactor ml — m" outside of the 
brackets we realize that this is a consequence of the quadratic dispersion relation connecting kinetic energy and wave 
number. This sum m' + m" leads eventually to the integer value n determining the steepness of the spacetime lines 
defined in Eq. (|22|). Hence many values of the summation indices m' and m" lead to the same steepness n _1 . 

According to Eq. ( pS] ) each J-function has a complex valued prefactor. This prefactor contains the term (— 1)"'' 
as well as the exponential exp[±i7rn(x' — x")/(2L)] or exp[±iwn(x' + x")/(2L)]. Again as in the argument of the 
(5-function either the sum or the difference of the integration variables x' and x" enters. Hence there is an interesting 
factorization property in K: When the difference x' — x" occurs in the argument of the ^-functions the phase of the 
complex amplitude contains the sum x 1 + x" and vice versa. This factorization property is essential for the emergence 



of the intermode traces as we discuss in Sec. [V 



We conclude this section by noting that the kernel K for the probability to find the particle at time t at the position 
x in the box consists of ^-functions aligned along straight spacetime lines. These world lines have discrete steepness 
and enter or leave the spacetime strip at x = at times determined by the sum and the difference of the integration 
variables x' and x" of the initial wave packet. 



IV. EMERGENCE OF INTERMODE TRACES 



In the preceding section we have realized that the kernel factorizes into sums of products of two functions. Each 
of the functions either depends on the sum x' + x" or the difference x' — x" of the integration variables x' and x", 
that is, of the position variables of the initial wave packet f*{x') and <p(x"). We are therefore lead to consider wave 
packets which satisfy the factorization property 

<p* (a/) <p{x") = <p+ (x' + x") <p- (x' - x") (28) 

where (p + and ip_ are new wave functions. Moreover, the initial wave packet has to satisfy the boundary conditions. 
Are there any wave functions that factorize and vanish at the walls? Here we do not want to go into a full discussion 
of this question but only mention that a narrow Gaussian satisfies both requirements in an approximated sense as 
discussed in Sec. [v|. 

For the case of factorizable wave functions we can now perform one of the two integrations with the help of the 
5- functions in the kernel K . Then the probability to find the particle at time t at position x reads 

W(x,t)= dx' dx" ip-{x' -x")<p+(x' + x")K{x,t\x',x"). (29) 



J J 

Under the assumption of a very narrow initial wave packet, it is possible to extend the integration limits to minus 
infinity and plus infinity. Motivated by the factorization property of the kernel discussed in the preceding section we 
now introduce the more convenient integration variables 

Z + = ±.(x' + x"), (30) 

and 

Z_ = -L(x'-x"). (31) 
We make use of the expression, Eq. (E3J), for the kernel K and the probability distribution becomes 
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W(x,t) = ~ ]T {-IT 1 de_/ ^ + ^(2L£_)^+(2L£ + ) 

x {e-<- 5 [ Xni i(x, t) - £+] + er™<- 5 [ Xn ,l(x, t) + £+] 

-e^S[ X nAx,t)-C-}-e- m ^d[xnA(x 1 t)+^}} . (32) 

The presence of ^-functions allows us to integrate over one variable. The remaining integral is the Fourier transform 

/• oc 



1 f°° 

0±(Kn) = ^J dye lhnV tp±(y) 
of the factorized wave packets and ip + to the wave number variable K n = nir/(2L). Hence wc find 

W(x,t) =- (-l) n ' i {^-(«n)¥'+[2LXn > K^*)]+0-(-'«r i )^+[-2iXr 1> i(^*)] 



(33) 



..n— — oo 



ip + (K n ) (p^[2Lx„j(x,t)} - tp + (-K n ) ip-[-2Lx„,i(x,t)}} 



When we change in the second and forth term the summation over n and I by defining n' = 
from the definition Eq. ( p2| ) of \n,i the relation 

X- n -l(x,t) = -XnA~x,t) 

we arrive at the equivalent formulation 

W(x,t) = - (-ir l {0-^n)^ + (2LxnAx,t)) + ^+(2L Xnd (-x,t))} 



and V 



(34) 
-I and note 

(35) 



1.71—- — OO 



- t P+( K n) [<P- (2L Xnil (x, t)) + </3_ (2LXn,l(-X, <))]} 



(36) 



This result is valid for any factorizable product ip*(x') ip(x") of initial wave packets. It brings out most clearly the 
way in which the initial conditions determine the visibility of the intermode traces. Indeed, according to Eq. (|36| ) the 
probability distribution W(x, t) to find the particle at time t and position x is a superposition of wave functions 
and ip+, evaluated along the spacetime lines Xn,i(x,t) and Xn,i{~ x ,t). The slope n _1 of these world lines is selected 
by the Fourier transforms (f>- and which represent the momentum distribution associated with the factorized 
wave packets. Hence, the factorized wave functions, </?_ and 9?+, and their Fourier transforms, 0_ and play the 
role of weight factors for the spacetime lines Xn,i(x,t) and Xn.i{~ x,t). 



V. EXAMPLE 



We now illustrate the result, Eq. (34), for the probability distribution W(x,t) by applying it to the special case 
of a Gaussian initial wave packet. Such a packet does indeed enjoy the factorization property Eq. (^8|). However an 
arbitrary Gaussian does not satisfy the boundary conditions enforced by the walls. For this reason we now focus on 
a Gaussian wave packet that is initially much narrower than the width of the box. 

We consider a Gaussian wave packet 



g(x) = (VttAx) 



-1/2 



exp 



1 1 x — x 

2 { Ax 



exp [ip(x — x)/h] 



(37) 



of width Aa; centered at x = x and moving with average momentum p = Tik > 0. Indeed this wave packet satisfies 
the factorization condition 



g*(x')g(x") = V+ (x'+x")ip-(x'-x") 
where the new factorized wave functions 



<P+(V) = (v^Aa 



,-1/2 



exp 



( y/2- 

\ Ax 



(38) 



(39) 



G 



and 



<P-(y) = (%/^ A ) 1/2 exp 



fJLV 

\2AxJ 



exp(-ipy/h) 



(40) 



are also Gaussian wave packets. 

After calculating the corresponding Fourier transforms and inserting them into Eq. ( |34| ) for the probability distri- 
bution W(x, t) to find the particle at time t and position x, one arrives at the following expression 



1 00 ( 

W(x,t) = — £ (-I)""' exp 



l,n=— oo 



■ exp 
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exp 




\Ak) _ 




\ Ax ) _ 



cos [2(kL\ n ,i - K n x)] > , 



(41) 



where Ak = 1/ Ax. 

We can identify in Eq. (|4l]) two different kinds of structures among all the possible intermode traces. 

We start our analysis from the Gaussian exp[— (k„ — k) 2 / Ak 2 } which is centered around the average wave number 
k. It takes on dominant values only for \n n — k\ ~ Ak. The other weight factor, exp[— (Lx n ,i — x) 2 / Ax 2 )], selects 
the spacetime lines starting at t = from the initial average position x of the particle in the box and from its images 
x + IL, periodically located along the x-axis. The combined effect of these two factors is to restrict the intermode 
traces to those connected with the classical motion. They therefore have a slope n _1 related to the average momentum 
and emerge from the initial location of the particle in in the potential well. 

A similar argument applies to the second term of Eq. ([li"|). In this case the spacetime lines Xn,i{x, t) are characterized 
by a negative slope — \ri\~ 1 satisfying the condition \n n + k\ ~ Ak. They originate from — x + IL. 

The third term of Eq. (|4l]) has the form of an interference term, weighted by two Gaussians. The Gaussian containing 
the quantized wave vector k„ is centered around zero. Therefore now the lines with small n, or equivalently with large 
slope are enhanced. According to the weight factor exp[— (Lx n ,i/ Ax) 2 ] these structures emerge at the time t = 
from the walls of the box, located at x = and x = L, and from their replicas periodically placed at x = IL. These 
are indeed the relevant intermode traces giving rise to the regular pattern in the spacetime plot of the probability 
distribution W(x,t). 

The last factor, cos[2(kLx n ,i ~ k„x)], provides a modulation depending on the initial conditions, that is average 
position and wave number. In the extreme case of zero momentum, shown in Fig. [l], the cosine factor reduces to 
cos(nirx/L). Hence whenever its argument is equal to ±7r/2mod(27r) the corresponding intermode traces are washed 
out. Indeed, in the example considered in Fig. [l the lines t/(T/2) = ±l/2(x/L — I) are not present because the cosine 
factor cos(n7r/4) vanishes for n — ±2. 




FIG. 1. Density plot of the probability distribution W(x,t) to find the particle at time t at position x. The initial wave 
packet is a Gaussian with width Ax = L/20 and average wave number k = centered at x = L/4. Darker areas correspond to 
minima of probability, while brighter ones represent maxima. 
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In the opposite limit of large average momentum, the cosine term is responsible for the fine structure of the quantum 
carpet. A comparison between Figs, jl] and ^ shows that by increasing the initial momentum, both canals and ridges 
develop neighboring structures: Consider, for example, the lines running along the main diagonals. In the case of 
non-vanishing initial wave numbers they split into an alternating series of canals and ridges. The detailed spatial 
behavior of such intermode traces is shown on top of Fig. @. 




FIG. 2. Density plot of the probability distribution of the particle in a box as a function of position and time. The initial 
wave function is a Gaussian of width Ax — L/20 and average wave number k = 20ir/L centered at x — L/4. On top we show a 
magnified view of the interference term in Eq. ([fl]) for n = 1, that is for the diagonal connecting the lower left corner and the 
upper right corner, and for the other main diagonal corresponding to n = — 1. Note that in the first case two ridges surround 
a canal, while in the second case the situation is reversed and two canals surround a ridge. 

VI. CONCLUSIONS 

We have studied the probability distribution of a particle in a one-dimensional box as a function of position and 
time. The corresponding density plot in spacetime is characterized by a regular pattern of intermode traces, that is 
canals or ridges cutting through and building onto an almost uniform background of probability. 

In order to explain this phenomenon we have related the probability distribution to the product of the initial wave 
packet and its complex conjugate via the kernel K. This formulation allows us to distinguish between the intrinsic 
characteristics of the system and the influence of the specific initial conditions. 

Indeed we have found that the kernel, which consists of the product of two Green's functions, contains already the 
spacetime structures. It can be expressed as a superposition of (5-functions whose argument is a linear function of 
time and position. The kernel is non-vanishing only along this set of infinitely many lines with discrete steepness. 
However, the initial conditions can either enhance or suppress the visibility of such structures. They therefore select 
the ones with the largest degeneracy from all possible intermode traces. 

We have shown that an initial wave packet which satisfies a certain factorization property is especially suitable 
to bring out the pattern in the probability distribution. An example of such wave functions is the Gaussian wave 
packet, discussed here. We emphasize however that the treatment presented here is more general. It even allows us to 
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investigate the case of a uniform initial distribution with discontinuities, analyzed by Berry in the context of fractal 
properties of the spacetime probability distribution. 
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APPENDIX: NEW REPRESENTATION OF V 

In this appendix we reexpress the function 

oo 

T>(r]; 0= Y exp[i7r(m' + m")rf\ exp {-iir{m' - m") [£ - (m! + m")2r + £]} (42) 

m' ,m" — — oo 

defining the kernel K, Eq. (|l9|), in terms of infinitely many 5- functions aligned along straight lines in spacetime. 
Throughout this appendix we use the dimensionless position variables £, 77 and £ normalized to the length of the box 
and the dimensionless time r scaled with respect to the revival time T. 

We note from the definition of T> that only the sum m! + m" and the difference m! — m" of the summation indices 
ml and m" occur. This suggests to introduce new summation indices m = m! + m" and k = m! — m" . However, this 
definition implies that m! — (k + m)/2 and m" = (m — k)/2. Hence when one of the two new indices is odd and the 
other is even the quantities m! and m" are not integer anymore. Therefore, m and k must both be either even or odd 
leading to the substitutions 

2m = m + m and 2k = m! — m" (43) 

or 

2m + 1 = m + m and 2k + 1 = m — m , (44) 

and to the representation 

00 

V= e n(2m)ri exp {-m{2k) [£_ - (2rn)2r + (}} 

k.77l — — OO 
OO 

+ Y e iTr(2m+1)l ' exp{-i7r[f-(2m+l)2T + C]} exp {-i7r(2fc) - (2m + l)2r + C]} . (45) 

k.rn— — 00 

When we recall the relation 

00 00 

Y exp(-2irik6) = S (° ~ ( 46 ) 

k— — 00 I — — 00 

we can perform the summation over k and arrive at 

00 

v= J2 e l < 2m ^ r ' 8[S, ^ {2m)2r + C, - l] 

/,m— — OO 
OO 

+ ^ (2m+1)ri exp {-iir[£ - (2m + l)2r + (}} 5[£ - (2m + l)2r + ( - 1} (47) 



l,m= — 00 

which with the help of 







simplifies to 



oo 

V= 2 e in(2m) vs ^_ (2m)2 T + (-l] 

l,m— — oo 

oo 

+ ^ e iT(2m+l), e -i,rl ^ _ ( 2m + 1)2t + C - I] 



(48) 



(49) 



l,m= — oo 

We can combine the sum over (2m) and the one over (2m + 1) into one sum over n when we recall the relation 

e~ M = (-1)' (50) 

and note that 



(-l) r 



1 for n = 2m 
for n = 2m + 1 . 



This yields the new representation 

oo 

T>(W,0= E (-irV™M[£-n(2r)+C-J] 



(51) 



(52) 



n,t— — oo 



which is central to the analysis of the kernel K in Sec. Ill 



[1] See for example D. Bohm, "Quantum theory" (Prentice-Hall, Englewood Cliffs, New York, 1951). 

[2] See for example G. Herzberg, "Electronic spectra and electronic structure of polyatomic molecules" (R. E. Krieger Pub. 
Co., Malabar, Fla., 1991). 

[3] See for example A. Bohr and B. R. Mottelson, "Nuclear structure" (Benjamin, Reading, Mass., 1969). 
[4] W. Kinzel, Phys. Bl. 51, 1190 (1995) and the letter to the editor by H. Genz and H.-H. Staudenmaier, Phys. Bl. 52, 192 
(1996). 

[5] M. V. Berry, J. Phys. A 29, 6617 (1996); M. V. Berry and S. Klein, J. Mod. Optics 43, 2139 (1996). 
[6] P. Stifter, C. Leichtle, W. P. Schleich, and J. Marklof, Z. Naturf. 52 a, 377 (1997). 

[7] J. Marklof, "Limit theorems for theta sums with applications in quantum mechanics" (Shaker Verlag, Aachen, 1997). 
[8] F. Grofimann, J.-M. Rost and W. P. Schleich, J. Phys. A 30, L277 (1997). 

[9] I. Marzoli, O. M. Friesch and W. P. Schleich, in Proceedings of the Fifth Wigner Symposium, ed. P. Kasperkovitz (World 
Scientific, Singapore, in press). 

[10] A. E. Kaplan, P. Stifter, K. A. H. van Leeuwen, W. E. Lamb, Jr. and W. P. Schleich, Physica Scripta, in press (1997). 
[11] D. L. Aronstein and C. R. Stroud, Jr., Phys. Rev. A 55, 4526 (1997). 

[12] P. Stifter, W. E. Lamb, Jr. and W. P. Schleich, in Proceedings of the Conference on Quantum Optics and Laser Physics, 
ed. L. Jin and Y. S. Zhu (World Scientific, Singapore, 1997). 



10 



